
clear
use "./gdp.dta"


gen time = quarterly(yq, "Yq")
tsset time
gen lngdp = ln(realgdp)
*gen dlngdp = D1.lngdp

tsfilter hp hplngdp = lnrealgdp
qnorm hplngdp //if year >= 1995 & year <= 2015, ///
mlabel(yq) mlabsize(vsmall) msize(vsmall) //

*& year!=2008 & year!=2009
*qnorm hplngdp
